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r S^ Abstract 

We study bipartite entanglement in systems of N identical bosons distributed in M 
different modes. For such systems, a definition of separability not related to any a 
priori Hilbert space tensor product structure is needed and can be given in terms of 
commuting subalgebras of observables. Using this generalized notion of separability, 
we classify the states for which partial transposition turns out to be a necessary and 
sufficient condition for entanglement detection. 

1 Introduction 

The characterization and quantification of non-classical correlations have become one of the 
most important tasks in quantum physicsQ Quantum correlations, besides allowing the 
implementation of classically unavailable protocols in information theory, are emerging as a 
fundamental tool in explaining the behaviour of many quantum systems, specifically those 
containing a large number of elementary constituents for which conventional approaches fail 
to provide a satisfactory description [21 [3] . 

Most of the studies devoted to the understanding and measure of the entanglement present 
in many-body systems have been focused on spin systems and the like, i.e. on systems where 
the microscopic constituents can be singly addressed. In this case, entanglement theory is 
well developed and various tests able to recognize and quantify non-classical correlations in 
such systems states have been identified and widely used in applications [U H] . 

The situation is less definite and clear when dealing with many-body systems made of 
identical particles, where the microscopic constituents can not be singly addressed and their 



1 See [T] and references therein. 
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properties measured: only collective, global operators are in fact admissible observables [5] [6]. 
Such systems are becoming more and more relevant in many applications, thanks to the 
recent advances in quantum optics, confined ultracold gases and superconducting systems, 
and more in general in the physics of mesoscopic systems and quantum devices. 

For systems in which the microscopic constituents are indistinguishable, the usual defini- 
tions of separability and of entanglement are inapplicable since the natural particle tensor 
product structure on which these notions are based is no longer available. More specifically, 
in a system made of N distinguishable constituents (e.g. qubits, spins, quantum dots...), 
a generic state u (i.e. a positive linear functional on the Banach space of the observables) 
is said to be entangled if it can not be written as a convex combinations of product states, 
namely 

co = ^Pk c*4 <8> <8> • • • <8> oj[ N) , Pk > , y^Pk = 1 , (1) 
k k 

where represents a state for the i-th constituent, with the associated "single-particle" 
Hilbert space / HS l \ This definition stems from the natural splitting of the total Hilbert space 
H of the system in the tensor product of the single constituent Hilbert spaces: 

u = n w ® n {2) ® . . . ® u {N) . (2) 

On the other hand, systems made of identical constituents can not be described by states as 
in (pQ) since the Hilbert space H must now contain only elements that are even (for bosons) 
or odd (for fermions) under any permutation of the elementary constituents. This fact 
points to the need of a generalized notion of separability and entanglement, not based on the 
tensor decomposition in (j2J), or more in general on the "particle" aspect of first quantization. 
The emphasis should shift from the set of system states to the corresponding algebra of 
observables and to the associated correlation functions, i.e. to a view closer to the second 
quantization description of many-body systems [7j-[T5]. 

Following this framework, a new definition of separability has been proposed in [T3j that 
applies equally well to all situations and reduces to the one given in flTJ) for systems of 
distinguishable particles. Although the notion of entanglement in many-body systems has 
been widely discussed in the recent literature [K>]-[27j, only limited part of those results are 
really applicable to the case of identical particles. We stress that the definition introduced 
in [13] is instead very general, suitable for all physical situations; in particular, it has been 
applied to the study of trapped ultracold bosons, leading to novel experimental testable 
predictions in quantum metrology [HI EE] . 

Having a general, consistent notion of entanglement is however not enough for a practical 
use in applications; as in the case of the standard definition based on the tensor decomposition 
(|21), one needs to identify criteria able to detect non-classical correlations through easily 
implementable tests [H [28] . 

In the case of bipartite entanglement in systems of distinguishable qubits, the operation of 
partial transposition has been identified as one of such criteria [291 [30]: the negativity of the 
partially transposed state is a sure sign of the presence of quantum correlations. Although 
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in general a sufficient condition for bipartite entanglement, the partial transposition criteria 
becomes also necessary in lower dimensions, either for a2x2-or 2x3- dimensional Hilbert 
space [29l [30] , or, in the case of continuous variable models, for bipartite Gaussian states in 
which one part contains just one mode [211 E2]- 

In the following, we shall discuss to what extent these results can be extended to bosonic 
systems made of a fixed number of elementary constituents. Using the generalized definition 
of separability and the language of second quantization, we shall classify the states for which 
the partial transposition operation results a necessary and sufficient criterion for detecting 
bipartite entanglement. As a byproduct of these results, a complete characterization of the 
states that remain positive under the operation of partial transposition will also be given, 
a fundamental step towards the full classification of bound entangled states in such bosonic 
systems. 

2 Entanglement and identical particles 

As remarked in the Introduction, for a generic many-body system made of distinguishable 
particles the usual notion of separability, as expressed by makes use of the particle 
tensor product structure ([2]) of the system Hilbert space, a clear heritage of the intrinsic 
addressability of the single elementary constituents: it reflects the natural multi-partition of 
the system into its elementary parts. 

On the other hand, states of the form ([I]) are not allowed states for systems made of 
identical particles; in fact, assigning a state a/ 1 ) to the first particle, co^ to the second and 
so on, would imply the possibility of distinguishing them through their states. Instead, the 
rules of quantum mechanics assign to a system of identical particles a total Hilbert space 
"H made of elements that are totally symmetric (for bosons) or antisymmetric (in the case 
of fermions) under any permutation of the particles; similarly, also the system observables 
must correspondingly be either even or odd under the same class of permutations [5j [6]. 

As a consequence, the notion of separability based on the decomposition ([T]) is no longer 
applicable in the case of bosonic (or fermionic) many-body systems: a suitable generalization 
encompassing it is needed. The basic idea is to give emphasis to the set of the system 
observables instead of that of its physical states; these two "points of view" are "dual" to 
each other: the connection between the two is given by the expectation value map that 
allows to express the average of an observable O as the value taken by the system state u on 
it, (O) = w((9)H Recalling [13], we shall first give an abstract definition of this generalized 
notion of separability and then apply it to the specific case of multi-mode bosonic systems. 

Let us consider a many-body system described by the Hilbert space "H and denote by 
i3("H) the algebra of all bounded operators on it; the observables of the system are part of 
this algebra. We shall introduce the notion of (bipartite) separability by considering couple 

2 A standard representation of this expectation value map is given by the trace operation over density 
matrices. 
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of commuting subalgebras of B(H) instead of focusing on partitions of its Hilbert space. We 
then introduce the following preliminary definitions: 

Definition 1. An algebraic bipartition of the algebra B(T-L) is any pair (Ai, A2) of commut- 
ing subalgebras of B(7i), A\,Ai C B(H); this means that any element of the subalgebra A\ 
commutes with each element of Ai, in short: [A\,A2] = 0. 

The two subalgebras Ai and A2 need not reproduce the whole algebra B(T-L), i.e. in general 
Ai U A2 C B{K). In this respect, the term "bipartition" used above should be interpreted 
in loose terms; nevertheless, in the cases discussed below, the considered mode partitions 
actually generate the whole algebra B{l-L). The notion of algebraic bipartition allows defining 
the system local observables: 

Definition 2. An element (operator) of B(H) is said to be local with respect to a given 
bipartition (A\, A2) if it is the product A\A 2 of an element Ai of A\ and another A 2 in „4 2 - 

We can now introduce the generalized notion of separability and entanglement: 

Definition 3. A state uj on the algebra £>("H) will be called separable with respect to the bi- 
partition (Ai, A2) if the expectation oj(AiA 2 ) of any local operator A\A 2 can be decomposed 
into a linear convex combination of products of expectations: 

u{A l A2) = Y,^\A l )ujf\A2) , A fe >0, ^A, = l, (3) 

k k 

where and are states on Bi^-C); otherwise the state u is said to be entangled with 
respect the bipartition {A\,A-i). 

This generalized definition of separability results meaningful in all situations and can be 
easily extended to the case of more than two partitions, by an appropriate, straightforward 
generalization of Definition 1 and Definition 2; specifically, in the case of an n-partition, 
Eq.Q would extend to: 

u(A 1 A2---A n ) = J2^^\A 1 )J*\A2)---J*\A n ) , A fe >0, ^> = 1 . (4) 

k k 

Further, when dealing with systems of distinguishable particles, Definition 3 gives the stan- 
dard notion of separability (cf. Eq.([T])). In this case, the partition of the system into its 
elementary constituents induces a natural tensor product decomposition not only of the 
Hilbert space 7-L as given in ([2]), but also of the algebra B(7i) of its operators; a direct appli- 
cation of the condition (J3J) to this natural tensor product multipartition immediately yields 
the decomposition ([T]). 

In this respect, it should be noticed that when dealing with systems of identical particles, 
there is no a priori given, natural partition, so that questions about entanglement and 
separability, non-locality and locality are meaningful only with reference to a specific choice 
of commuting algebraic sets of observables; this general observation, often overlooked in the 
literature, is at the basis of the definitions given in ([3]) and (T4]). 
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A special situations is represented by pure states. In fact, when dealing with pure states 
instead of general statistical mixtures and bipartitions that involve the whole algebra B(T-L), 
the separability condition in ([3]) (and similarly for simplify, becoming: 

oj(A 1 A 2 )=J 1 \A 1 )J 2 \A 2 ) ; (5) 

in other terms, separable, pure states are just product states. Indeed, recall that a state uj is 
said to be pure if it can not be expressed as a convex combination of other states; therefore, 
the sum in must contain just one term, otherwise it would be in contradiction with the 
hypothesis of uj being pure. Further, by taking either A\ or A 2 to be the identity operator, 
from the expression in (JS} one sees that the factor states uj^ and uj^ coincide with the 
restrictions of uj to the first and second partition, respectively, so that ([5]) reduces to: 

uj(A x A 2 ) =u(Ai)u){A 2 ) ■ (6) 
This observation will be useful in the following. 



3 Multimode boson systems 

Henceforth, we shall focus on bosonic many-body systems, whose elementary constituents 
can be found in M different states or modes. This is a very general framework, useful for the 
description of physical systems in quantum optics, in atom and condensed matter physics. 
Among the many models fitting this paradigm, those consisting of ultracold gases of bosonic 
atoms confined in multi-site optical lattices is of special relevance. They turn out to be a 
unique laboratory, both theoretically and experimentally, for the study of quantum effects in 
many-body physics, e.g. in quantum phase transition and matter interference phenomena, 
and also for applications in quantum information (e.g. see [33]- [H3], and references therein). 

In dealing with (bosonic) many-body systems, it is most appropriate to adopt the second 
quantized formalism. Let us thus introduce creation a] and annihilation operators a,, i = 
1,2,..., M, for the M different modes that the bosons can occupy; they obey the standard 
canonical commutation relations, [a,, a]] = and generate out of the vacuum state |0) 
single-particle orthonormal basis states \i) = a||0). The total Hilbert space H of the system 
is then spanned by the many-body Fock states, 

\ni,n 2 , ...,n M )= ] (a})™ (4)- . . . (a^f - |0> , (7) 

where the integers rii,n 2 , . . . , tim are the occupation numbers of the different modes. In this 
language, symmetrization of the elements of "H, as required by identity of the particles filling 
the M modes, is automatically guaranteed by the commutativity of the creation operators. 

Furthermore, notice that the set of polynomials in all creation and annihilation operators, 
{aj, di}, i = 1,2,...,M, form an algebra that, together with its norm-closure, coincides 
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with the algebra B(H) of bounded operators of the M-mode bosonic systemjl 

A bipartition of this algebra is determined by an integer m, < m < M, splitting the 
collection of creation and annihilation operators into two disjoint sets {a\, a* \i — 1, 2 . . . , m} 
and {a^, a a , \ a = m + 1, m + 2, . . . , M}; all polynomials in the first set (together with 
their norm-closures) form a subalgebra Ax, while the remaining set analogously generates 
a subalgebra A%. Since operators pertaining to different modes commute, one sees that 
[Ai, Aq\ = and thus the pair (Ax, Az) indeed forms an algebraic bipartition of the whole 
algebra For sake of clarity, we shall use latin indices to label operators in the set Ax, 

while adopting greek indices for those in A2 and similarly for the mode occupation numbers 
(see below). 

As explained in the previous Section, from an algebraic bipartition (Ax, A2) one inherits 
the notion of locality for the elements of B(Ji) (see Definition 2) and consequently that of 
separability and entanglement (see Definition 3). 

The case in which the two commuting algebras Ax and A2 are generated only by a subset 
M' < M of modes can be similarly treated. Indeed, all operators in B(TL) pertaining to the 
modes not involved in the bipartition commute with any element of the two subalgebras Ax 
and A2, and therefore effectively act as "spectators" . As a consequence, all the considerations 
and results discussed below holds also in this situation, provided one replaces the total 
number of modes M with M', the actual number of modes used in the chosen bipartition. 

Given our bosonic, M-mode many-body system, we shall first consider pure states and 
give a complete characterization of the separable ones, through the following: 

Proposition 1. A pure state is separable with respect of the above (Ax, A2)- 
bipartition (i.e. it is (Ax, A2) -separable) if and only if it is generated out of the vacuum 
state by a (Ax, A2)-local operator, in short 

\^)=V(al...,al)-Q(al +1 ,...,a{ I ) |0> , (8) 

where V, Q are polynomials in the creation operators relative to the first m modes, the last 
M — m modes, respectively. By Definition 3, a pure state that can not be cast in the 
form (0) is therefore (Ax, A2) -entangled. 

Proof: First of all recall that in the present situation the condition of separability reduces to 
the simpler expression clearly, the state in (|S]) satisfies it by taking for the expectation 
value of a generic operator A G B(%) the usual state-average: u(A) = (i/;\A\i/;). In order to 
prove the converse, i.e. that from the separability condition (JS]) the expression (jS]) follows, 

3 Strictly speaking, the polynomials themselves are not bounded operators (nor are the single creation 
and annihilation operators); the standard way to properly define them is through the introduction of the 
Weyl operators and the corresponding algebra of bounded operators, from which polynomial operators can 
be defined by differentiation [2U IH] • 

There is no loss of generality in assuming the modes forming the subalgebras Ai (and A2) to be 
contiguous; if in the chosen bipartition this is not the case, thanks to the commutativity of operators 
pertaining to different modes, one can always relabel the modes in such a way to achieve such convenient 
ordering. 
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let us decompose \ip) in the Fock basis given in (J7|): 

\4>) = l ra i' ra 2, • • • ,n M ) , 5^|C{ n }| 2 = 1 , (9) 

{n} {n} 

or equivalently, taking into account the bipartition (m, M — m) of the modes, 

W) = ^ C {k}>{a y\k 1 ,...,k m )\a m+1 ,...,a M ), (10) 
«,{<*} 

where the sums are over all integers hi, i = 1, . . . , m and aj, j — m + 1, . . . , M. Insert- 
ing this decomposition in the separability condition ([6]) and choosing A\ = (aj)* 5 ' 1 • • • (aj„) fe " 1 • 
(ai) kl ■ ■ ■ (a m ) km and A 2 = (a m+ i) a ' m+1 ■ ■ ■ (aj u ) a ' M -(a m+ i) a ™ +1 • • • (a M ) aM , one gets the follow- 
ing set of constraints on the complex coefficients C (the bar signifies complex conjugation): 

C{k'},{a'} C{k},{a} = C {k},{P} ^2 C W>W} G {l},{a} " i 11 ) 

Taking k' = k and a' = a, this expression becomes 
|C{fc},{a}| 2 = D {k} D' {a} , D {k} = ^ \ C {k},{P}\ 2 ^ > D {a} = E I^W.wf - ' ( 12 ) 

{1} 

so that the modulus of the coefficient C factorizes, and one can write: 

%},W = (D {k} ) 1/2 (D' {a} ) 1/2 e»«.M , D W = 1 > E D U = 1 • ( 13 ) 

{k} {a} 

Insertion of this expression back into the condition (fTT]) produces additional constraints for 
the phases 6: 

e i(<?{fc},M-0 {fe / }:{Q / } ) _ e *( e {fe},{^}- 6, {fc'},{/3}) e i ( e {i} i {Q}- 6, {;},{«'>) . (14) 

From this expression one immediately sees that the phase difference (9^},{i3} — @{k'},{/3}) 
does not depend on the labels {/?} and vanishes for k = k', while (0{i},{ a } — #{i},{ a '}) is 
independent from the labels {1} and vanishes for a = a'. With this information, one can 
rewrite the decomposition ffTUj) as 

iv) = j2( D w) 1/2 \ k ^---i k ^ E(%) 1/2e%Wo} i a ^'->^> ( 15 ) 

{fc} {a} 

{k} {a} 

where {a'} is any given set of reference occupation numbers. However, as observed before, 
the difference (0{k},{ a } — 6{k},{a>}) is independent from {k}, so that the expression of \ip) given 
above is precisely of the form ([8]). □ 
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When the state of the bosonic many-body system is not pure, it can be described by a 
density matrix p, that, using the Fock basis of (JZ]), can be written as: 

P= ^2 P{n},{n'} \ni,n 2 ,...,n M )(n' 1 ,n' 2 ,...,n' M \ , ^p {n}i{ „ } = 1 . (16) 

{n},{n>} {n} 

The set of all states form a convex set, whose extremals are given by the pure ones. Propo- 
sition 1 can then be used to characterize separable mixed states: 

Corollary 1. A mixed state p as in / fig)) is {A\, A2) -separable if and only if it is the 
convex combination of projectors on pure (Ai, A2) -separable states; otherwise, the state p 
is (Ai, A2) -entangled. 

Until now, we have considered many-body systems formed by an arbitrary number of 
elementary constituents; we shall now restrict the discussion to systems with a given number 
iV of bosonic particles. In this case, the Hilbert space H becomes finite dimensional: it is 
still spanned by the Fock states ((7J); however, since the total number operator Y^=i a ! a « * s a 
conserved quantity, the occupation numbers must satisfy the additional constraint Y2i=i n i = 
N, i. e. all states need to contain exactly N particles. 

Note that in absence of this number "superselection" rule, the many-body Hilbert space 
H spanned by the Fock states fl7]) possesses a mode tensor product structure similar to the 
particle tensor product in (j2J), % = V,^ <8> Ti.^ ® ... <g) T-6- M \ where T-&' now represents 
the z-th mode Fock space spanned by the states \rii) ~ (a|) ni |0), rii G N. This structure is 
clearly lost when the number of the system constituents is fixed, since now all states must 
contain the same number of bosons. 

The conservation of the total particle number has consequences also on the form taken by 
pure separable states \ip); they still have to be in product form as in flH]), but now V and Q 
need to be homogeneous polynomials of complementary degree: if V is of degree k, than Q 
must be of degree N — k. In other terms, given the mode partition (m, M — m) leading to 
the bipartition (A±, A2) of the whole algebra B(T-L), if there are k particles filling the first 
m modes, then there must be N — k in the remaining M — m modes; only in this way the 
state \ip) will contain the prescribed N bosons. 

Examples of N bosons pure separable states are the Fock states; indeed, using the notation 
and specifications introduced before, they can be recast in the form (jSJ): 

m M 

\k 1 ,...,k m )\a m+1 ,...,a M ), J ^k i = k, ^ atj = N - k , < k < N ; (17) 

i=l j=m+l 

by varying the integers ki, aj and k, these states generate the whole Hilbert space 7i. This 
basis states can be relabeled in a different, more convenient way as: 

»r , a {k + m-l\ , fN-k + . 

\k,a;N-k,a'}, a = l,2,...,( fe J , a' = 1, 2, . . . , f ^ 
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as before the integer k represents the number of particles found in the first m modes, while 
a counts the different ways in which those particles can fill those modes; similarly, a' labels 
the ways in which the remaining N — k particles can occupy the other M — m modes. 
Clearly, due to the orthonormality of the states in f[T7|) . in this new labelling one has: 
(k, a; N — k, a'\l, r;N — I, r') = 5ki 5 aT 5 a i r >- Using this notation, a generic mixed state p can 
then be written as: 

N N 

P=^2 Pkaa',irr' \k,a;N - k,a')(l,r;N - l,r'\ , y] Pkaa'Ma' = 1 ■ (19) 

k,l=0 <t,<j',t,t' k=0 cr,cr' 

Different orthonormal basis of separable pure states can be chosen for % and will be used in 
the following; they are all characterized by the integer k that fixes the number of particles 
filling the first m modes, plus some additional indices able to distinguish orthogonal states 
with the same number of particles in the two partitions. 



4 The partial transposition criterion 

In the previous Section, we have seen that, given a bipartition (A±, A2) of the iV-boson 
algebra, obtained through the (m, M — m)-partition of the M modes, separable pure states 
containing k particles in the first modes (and N — k in second ones) are in product form as 
given in ([8]), or equivalently, when expressed in the basis of the Fock states (fT8l) : 

\^) = J2ci k) Ci k , ) \k,a;N-k,a') , (20) 

a,cr' 

while separable mixed states are convex combinations of their projectors. However, to de- 
termine whether a given density matrix can be written in separable form is in general a hard 
task. In the case of systems of distinguishable particles, easily implementable separability 
tests have been constructed. In the case of bipartite systems, one of the most useful criteria 
involves the operation of partial transposition [2"9l [SO] : a state p for which the partially 
transposed density matrix p is no longer positive is surely entangled. This lack of positivity 
can be witnessed by the so-called negativity [HH HU [1] : 

M(j>)= ^(llp||i-Tr[p]) , ||p||i=Tr 

Indeed, if p is not positive, then || p \\i> Tr[p] = Tr[p] = 1 and as consequence N(p) > 0. 
Unfortunately, this criterion is not exhaustive, since there are entangled states that remain 
positive under partial transposition. The negativity is nevertheless an exhaustive entangle- 
ment witness for systems composed by two qubits or one qubit and one qutrit [29j [30], or, 
in continuous variable systems, for two- mode Gaussian states [31]. 

We shall now examine to what extent these results, valid in the case of bipartite systems of 
distinguishable particles, can be extended to systems composed by N bosons; we will see that 
the identity of the particles induces novel features in the theory of many-body entanglement. 
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We shall work within the generic (A±, ^.2)-bipartition introduced above, induced by the 
(m, M — m)-splitting of the M modes. First notice that the operation of partial transposition 
{e.g. with respect to the first subalgebra) maps separable states into separable states. This 
is immediate for pure states as in fl2"01 . since in this case partial transposition amounts to 

the complex conjugation of the first coefficient, | — > C a ; being convex combinations 
of projection on pure separable states, also separable mixed states are similarly mapped 
into perfectly admissible density matrices. As a consequence, also in the case of bosonic 
systems, negativity of a partially transposed density matrix would definitely signal (Ai, A2)- 
entanglement. 

For a generic density matrix p as in ( 119|L the computation of the negativity N(p) involves 
the evaluation of the trace norm || p \\% of the partially transposed matrix: 

N 

P=J2 Yl Pk^'MT' \l,r;N-k, a') (k, a;N-l,r'\, (22) 

k,l=0 o-,ct',t,t' 

usually a rather difficult task (however, see Section 5 below). Nevertheless, there are relevant, 
quite general instances in which the evaluation of M{p) can be performed analytically. 

As a prelude to more complex situations, we shall first discuss the case of a system of iV 
bosons that can occupy just two modes, i.e. M = 2, each of the two algebras A± and A2 
being generated by the creation and annihilation operators in just one of the two modes; 
this example is quite important in applications since it is the paradigm for the description 
of ultracold bosonic gases confined in double well potentials [45, 46J. One finds that in this 
case the negativity is an exhaustive entanglement witness, as stated by the following: 

Lemma 1. A two-mode state of an N -boson system is entangled with respect to the 
bipartition into single modes if and only if its partial transposition is negative. 

Proof. A generic state for a two mode iV-boson system can be written as in (1191) . although 
in a simplified form: 

N N 

P = \k;N -k)(l;N-l\ , 53 p fc ,fc = l, (23) 

k,l=0 k,l=0 

since in this case a single integer is sufficient to label the iV + 1-dimensional basis of Fock 
states. The partial transposition operator (with respect to the first mode) maps this state 
to an operator p on a larger Hilbert space than that with a fixed N, 

N 

p=J2 PKi \l;N-k)(k;N-l\ , (24) 

k,l=0 

such that the combination 

N 

P ] P= Yl l^l 2 \k;N-l)(k;N-l\ (25) 

k,l=0 
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is diagonal. Its trace norm can then be readily computed, as the negativity in (121]): 



M(P) = g E M ■ ( 26 ) 

k^l=0 

This expression vanishes if and only if the state is diagonal in the Fock basis, but this is 
equivalent to separability, since in this case p results a convex combinations of projectors 
over Fock states. □ 
These conclusions can be extended to the case of bosonic systems in M modes by choosing 
a (Ai, ^.2)-bipartition for which the algebra Ai is generated by creation and annihilation 
operators of one single mode, while the remaining M — 1 modes generates A%. Also in this 
case the negativity results an exhaustive entanglement criterion as stated in: 

Proposition 2. A M-mode state of an N -boson system is entangled with respect to the 
the bipartition into one mode and the rest if and only if its partial transposition is negative. 



Proof. In this case, a generic state of the system can be written as in (Il9p. dropping 
however the first greek labels in the Fock states: 



N N 



P = E E Pka ' h \k', N ~ k,cr)(i,N - l,r\ , EE Pka ' ka = 1 " ( 27 ) 



k,l=0 a,T k=0 a 

The operation of partial transposition with respect to the first algebra can be readily imple- 
mented 

N 

P = E E \l;N-k,a) (k; N-1,t\ , (28) 

k,l = <T,T 

so that 

N 

?p= E E hk^r \^i)(r kl \ , (29) 

k,l=0 a,T 

where 

/ \ -1/2 

= El^l E PW \k;N-l,a') . (30) 

t ' a' 

By separating diagonal and off-diagonal terms in k and I, the positive matrix in (|29|) can be 
further decomposed as 

p^p = D + R, (31) 

d = E E (E W°m\ 2 ) \rkk)(<f>i k \ , (32) 

k=0 cr t ' 

R = E E (El^.r) 1^x^,1 • (33) 
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Thanks to the orthogonality of the vectors \<p%i) f° r different indices k or /, one finds that 
the products DR and RD vanish, and thus one has: 

7^5 = VD + VR. (34) 

In order to explicitly compute \/~D, we make a change of basis in the second M — 1 partition 
passing from the Fock states to another set of separable states adapted to p, such that 
its components along this new basis satisfy: p k a,kr = Pka,ka S aT ; notice that this is always 
possible through suitable local, unitary transformations diagonalizing for each k the matrices 
Ai&r = [/0fc<7,fer] • I n this new basis, the vectors \<j>%i) become orthonormal in all three indices, 
and as a consequence, recalling (I3"2"j) . 

N 

v^ = EE \p*m\ I I ; (35) 

k=0 (J 

further, since the diagonal elements Pkakc are non-negative, one has: Tr[v^D] = 
J2k=o ^2(7 PkaM = 1; being p normalized; it follows that the negativity (l2Tj) of the gen- 
eral state (1271) is simply given by 

Af(p) = ^Tr [v^R] . (36) 

Recalling the explicit expression of R in fl33|) . one immediately sees that the negativity 
vanishes if and only if the elements p kcr .i T for k ^ I are all zero. But this means that 
p, expressed in the new adapted basis, is the convex sum of projectors on separable pure 
states. □ 

In the above proof we have used the partial transposition operation with respect to the 
first party. The same conclusion nevertheless holds by choosing the partial transposition 
operation with respect to the second party, since it can be obtained as the combination of 
the partial transposition with respect to the first party followed by the total transposition, 
and this second operation does not change the negativity. 

The previous proof suggests the existence of a larger class of bipartite states for which the 
negativity of the partial transposed states uniquely signals the presence of entanglement. 

Corollary 2. Let us consider the class of (m, M — m) mode bipartite states / fTPj) whose 
principal minors p k(TCT i kTT i for each k turn out to be diagonal in a basis of separable pure 
vectors. Then, these states result entangled with respect to the fixed bipartition if and only 
if their partial transposition is negative. 

Proof. This statement can be proved through a straightforward generalization of the steps 
followed in proving Proposition 2. In the chosen basis, one can decompose a generic state p 
as in 0191) . Then, by direct computation, one arrives at a decomposition p^p = D + R similar 
to the one in (l3~H - fl33"l) : the class of considered states are such that p karj i kTT i oc 5 UT 5 a i T i, and 
this allows to conclude that Tr[\AD] = 1. Then, as before, one finds that the negativity 
M{p) = Tr[V r R]/2 is zero if and only if the off-diagonal elements pw,irT'; for k ^ I, vanish, 
i.e. if and only if the state p is separable. □ 
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5 States remaining positive under partial transposition 



Although the results proven so far cover a large set of bosonic states, the negativity of a 
partially transposed state is not in general a necessary condition for entanglement: there are 
states that remain positive under partial transposition (they are usually called PPT states), 
but nevertheless turn out to be entangled. As a first step towards their classification, it is 
important to give an explicit characterization of iV-boson states that are PPT; this can be 
obtained by means of techniques similar to those used above. 

Proposition 3. A general (m, M — m)-mode bipartite N-boson state l[W\) is positive 
under partial transposition (PPT) if and only if it is "block diagonal" , namely of the form 

N 

P = ^2pk , (37) 

fc=0 

with 

Pk= ^2 Pkaa',krr' \k,cr;N - k,a')(k,r;N - k,r'\ , (38) 

and its principal minors pk are PPT. 

Proof Given the decomposition ([19]) of a generic state p in the Fock basis adapted to the 
chosen (m,M — m)-bipartition and its partial transposition (I22j) . one sees that the positive 
operator p^p can be cast in the form: 



E E (El^-'f) i^ t >(^ t i , (39) 

t. 7 — rr -r \ rrl^l ' 



k,l=0 cr,r v ct't 

where the vectors 

/ \ -1/2 

l*«> = I EKwl ) E *Wr' \k,a';N-l,r') , (40) 

are orthogonal with respect to the latin indices, but in general not in the greek ones. Nev- 
ertheless, this is sufficient for proving that the following decomposition: 

N 

^ = J2^k + Vn, (4i) 

fc=0 

similar to the one encountered before in a simpler setting, holds; explicitly, one has: 

v * = E (E W'rMr'f) , (42) 

ct,t ^ ct't' ' 

n= E E (EK'^'I 2 ) hhdwi. ( 43 ) 

k^l = CT,T ^ CT'T' ' 
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so that one immediately checks that D^V^i = for k ^ k' and also T>k1Z = = lZT>k- 
Further, by direct computation, one finds that T>k = (pkYpk, where pk is the following 
operator built with the "block-diagonal" components of p (cf. fl38l) ): 



N 



Pk 



^ pkaa',krr' \k,a;N - k,a')(k,r;N - k,r'\ , ^Tr[p fc ] = 1 



k=0 



the second expression coming from the normalization of p. As a consequence, the negativity 
of the generic iV-boson state p can be written as: 

N 1 r i 

Af(p) = J2 Af (P*) + 2 TT V M \ ' (44) 

k=0 

in terms of the negativities Af(pk) of the operators p^. Since TZ as given in fj4^|) is the (direct) 
sum of positive operators, M{p) in (J33J) vanishes if and only if each term does, namely p is 
precisely of the form (1571) . and further M(pk) =0. □ 
This complete characterization of the iV-boson PPT states clearly simplify the search for 
bound entangled states; indeed, the identification of entangled bosonic many-body states 
that remain positive under partial transposition, in general a very difficult task, is reduced 
by Proposition 3 to the more manageable problem of identifying quantum correlations in 
their diagonal blocks (|38|) . for which standard techniques can be used pQ. 



6 Applications and examples 

The result in (|44p gives the expression for the negativity of a generic iV-boson state p, 
once written in terms of the Fock basis adapted to the chosen bipartition as in (fl9l) ; from 
its explicit form one can appreciate that the statements of Proposition 2 and Corollary 2 
proven before require much more stringent conditions than the ones assumed in the previous 
Proposition 3, specifically, that the negativity of all diagonal blocks of p be identically zero, 
M{pk) = 0. In this respect, the conclusions of Proposition 3 are therefore more general, but 
weaker, than the ones dealing with the partial transposition criterion in Section 4. 

Yet, the general expression (J33J) of the negativity gives some additional information on 
entangled iV-boson states in general; they can be summarized by the following: 

Corollary 3. If a generic (m, M — m)-mode bipartite state ( fTPj) can not be cast in block 
diagonal form as in l[3?\) . ( T38j) . or if it can, at least one of its diagonal blocks is not PPT, then 
the state is non-separable and its entanglement can be detected by the partial transposition 
criterion. 

In other terms, a given many-body systems is very likely to be found in a state having 
non-vanishing off-diagonal elements and therefore in an entangled state, since in this case 

Tr \pR, 7^ 0, giving rise to a nonvanishing negativity, M{p) > 0. Indeed, this is the 
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case in many physical situations involving quantum phase transitions in condensed and 
superconducting systems [UJ, where the so-called "off-diagonal long range order" states 
play a fundamental role (e.g. see [4~8]-[52]). 

Nevertheless, also block-diagonal states of the form (1571) . (13"5|) can be encountered in the 
study of many-body systems, due to the presence of noise. Indeed, actual many-body sys- 
tems, albeit carefully screened, can not be considered totally isolated from the external 
environment; its presence unavoidably produces noise and dissipation^ The most basic (and 
ubiquitous) effect induced by the environment is the so-called "dephasing noise", leading 
to decoherence, i.e. to the suppression of the off-diagonal elements in the system density 
matrix. 

More specifically, in our iV-boson setting, let us consider a general (m, M — m)-mode 
bipartition of our system, initially prepared in a generic state p: it can be decomposed as in 
flTTJ]) in terms of the Fock states \k, a; N — k, a') containing k bosons in the first partition and 
N — k in the remaining one. Using the mode creation and annihilation operators introduced 
at the beginning of Section 3, one can now consider the following collective observable 



M 



v = 22 a \ a i - ^2 a « aa ' ( 45 ) 



c*=m+l 



that counts the relative number of bosons in the two partitions; it clearly commutes with the 
total number operator J2iLi a l a i = N. The effects of dephasing noise can then be described 
by the following time-evolution master equation [53]-[57j: 



d P (t) 
dt 



l(v P (t)V- l -{v\p(t))^ , (46) 



where the positive coupling constant 7 measures the strength of the noise. Since the Fock 
states are eigenstates of V, V \k,a; N — k, a') = (2k — N) \k,a; N — k, a'), one easily finds 
that the matrix elements pkaa' ,i, T T'(t) = (k, a; N — k, a'\p(t)\l, r;N — I, r') obey: 

d 

Pkaa',l,rr'(t) = -7 { k ~ if Pkaa',l,rr'(t) ■ (47) 

It thus follows that this kind of noisy, irreversible dynamics indeed leads to the suppression 
of the entries that are off-diagonal in the partitions occupation number: 

N 

=H H e-^ fc -') 2 pw.Itt* \k,a;N-k,a')(l,r;N-l,r'\ . (48) 

k,l=0 a,a' ,t,t' 

As a consequence, in the long time regime, p(t) becomes diagonal with respect to this 
occupation number, the faster, the larger the coupling 7 is; at the end, p takes precisely the 
block-diagonal form given in (13T1) . (138]) . 



3 For an introduction to the theory of open quantum systems, see |53j-[57] 
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An experimentally relevant instance in which such diagonalization may rapidly occur is 
that of a system of N bosons with an internal (two-dimensional) hyperfine degree of freedom, 
confined in a double well potential]^ any boson can be found in one of the two wells with one 
of the two possible internal "polarizations" , so that in this case M = 4. A natural bipartition 
is the spatial one defined by the two wells (giving m — 2). The associated Fock basis states 
\k, a; N — k, a') are labelled by the integer k giving the number of bosons in the first well 
(N — k are in the remaining one) and by the additional integers a, a', e.g. labelling the 
number of bosons in the hyperfine ground state in the first, respectively second well; recalling 
([18]) . at fixed k, these states generate an Hilbert space isomorphic to C k+1 <g> C N ~ k+1 . 

The dephasing noise acts now on bosons in different wells, destroying the coherence be- 
tween the "spatial" modes, leaving instead untouched the internal, hyperfine ones. The 
resulting density matrix is block-diagonal, p = p k , where each block pk is a (non normal- 
ized) density matrix on C fc+1 ® C +1 . Clearly, as stated by Corollary 3, a non vanishing 
negativity of one of the p k s would signal entanglement for the iV-body state described by 
the total density matrix p. Furthermore, note that, thanks to Proposition 3, PPT states p 
that are entangled are also possible, provided N > 4: indeed, recalling the results of [29| 130]. 
only in this case the matrix p can accommodate diagonal blocks pk large enough to allow 
PPT entanglement. 

As an explicit example, consider the system composed by four bosons (N = 4), that can 
occupy four different modes (M = 4), and choose m = 2, so that each partition is made 
of two modes. The block diagonal density matrix p in ( 157)1 contains now five terms. Each 
diagonal block p k is fully described by the corresponding coefficient matrix p k(Ta , kTT r. for 
k — 0, 4 these matrices are 5x5,8x8 when k — 1, 3, while in the case k — 2 it turns out 
to be 9 x 9. Let us focus on this last p k =2 block: 

Pk=2= ^ Pvv'.tt' W;ct')(t;t'\ , (49) 

where on the r.h.s. we have suppressed for simplicity the labels k = N — k = 2. Since now 
cr, a', r, t' = 1, 2, 3, the Hilbert subspace spanned by the set of six Fock vectors {\a; a'), a, a' = 
1, 2, 3} is isomorphic to C 3 (g) C 3 , i.e. to a two qutrit space. For such a system, a bound en- 
tangled state has been presented in [59] and its 9 x 9 density matrix explicitly written in the 
C 3 (g) C 3 standard basis. By identifying the six Fock vectors {|cx; a')} with such a basis, one 
can choose for the coefficient matrix p acr > :TT i precisely that PPT entangled density matrix. 
Assuming that the remaining blocks p k , k = 1,3, are separable (the other two blocks p k , 
k = 0,4 are always separable), the density operator p = ^2 k=Q Pk results a bound entangled 
state for the system. 

As a final remark, notice that all above results can be generalized to the case of systems 
where the total number of particles is not fixed, but commutes with all physical observables 
(i.e. we are in presence of a superselection rule [61]). In such a situation, general density 

6 For actual realizations, see [551 IBU] . 
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matrices are incoherent mixtures of states of the form (1191) with fixed number N of particles: 

P = ^2 ^ NpN ^n>0, = 1 . (50) 

N N 

In this case, a suitable entanglement measure is the weighted sum of the negativities of each 
component p^, i.e. J2n ^nN(pn)- All the arguments discussed in the previous sections 
apply to the single pn, and therefore the statements of Lemma 1 and Propositions 2, 3 also 
hold in this more general settmgfj At the same time, the effects of the dephasing noise given 
by the dynamics fH6|) are still the asymptotic block diagonalization of each p^. 

Much more difficult is the treatment of the case of many-body states with fluctuating 
number of particles, i. e. of density matrices that are coherent mixtures of states with different 
N. At least when the number fluctuation is small, one can nevertheless advance some general 
considerations, based on the property of continuity satisfied by the negativity in fl2T]) . If one 
starts with a iV-boson state with strictly positive negativity (i.e. an entangled state), by 
the permanence of sign, perturbations by states with small number fluctuation should not 
be able to turn it into a separable state. On the contrary, any similar small perturbation 
of a separable iV-boson state (i.e. with vanishing negativity) will in general turn it into an 
entangled state, although with infinitesimal negativity; this is so because after perturbation 
the state will no longer be in block diagonal form. 

7 Outlook 

In the case of bosonic many-body systems composed by a fixed number of identical particles 
the usual notions of separability and entanglement are inapplicable since the Hilbert space 
particle tensor product structure on which they are based is no longer available. A generalized 
definition of separability valid in all physical situations can however be formulated: its based 
on commuting algebras of observables instead of the particle tensor decomposition of states. 
This new framework makes it apparent that in systems of identical particles, the notions 
of separability and entanglement are meaningful only in reference to a specific choice of 
commuting sets of observables. 

Using this generalized definition, we have studied bipartite entanglement in systems com- 
posed by N identical bosons that can occupy M different modes. More specifically, we have 
analyzed to what extent the partial transposition operation can be used to detect quantum 
correlations in such systems and found that in general it is a much more exhaustive criterion 

7 Being in presence of the number superselection rule, admissible pure states of the system need to be 
eigenvectors of the total number operator, while mixed states must be decomposed as their convex combina- 
tions. In general, the partial transposition operation mixes sectors with different N, giving rise to positive 
partially transposed matrices for states p not satisfying the separability criterion ([3]); indeed, this happens 
for states p whose separable convex decompositions always involve superpositions of different total number 
eigenspaces. Therefore, in presence of incoherent fluctuations of the number operator, the negativity as 
defined in (j!?Tj) docs not appear to be a useful entanglement measure. 
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for detecting bipartite entanglement than in the case of systems of distinguishable particles. 
Indeed, by giving explicit expressions to the negativity, the quantity that measures the lack 
of positivity of the partially transposed density matrix, we have shown that the operation of 
partial transposition gives a necessary and sufficient test of entanglement detection for very 
general classes of bosonic states. 

Even more relevant is the fact that such classes of states find direct applications in many 
physical models. We have already remarked that the two-mode case (M = 2) represents the 
paradigm for discussing quantum effects in the behaviour of ultracold gases of bosonic atoms 
confined in double-well traps; these systems have recently attracted a lot of attention in view 
of their metrological applications [6 2] - [65] , and in particular in relation to the possibility 
of constructing interferometric devices able to overcome the so-called shot-noise limit in 
parameter estimation [66]-[77]. 

Similarly, systems with several bosonic modes are used to model ultracold gases in optical 
lattices and other many-body condensed systems [2]- [I], [33]- [ID]) in particular in reference 
to quantum phase transitions [4T] . For such systems, the classes of entangled states discussed 
in Corollary 2, 3 find application as "off-diagonal long range order" states |48]-|52]. 

Even the block diagonal states analyzed in Proposition 3 are of relevance in actual real- 
izations. Indeed, as discussed in the last Section, they are the natural outcome of dephasing 
noise; this makes quite realistic the possibility of producing PPT entangled states in actual 
experimental apparata, e.g. those described in [581 ED]) through the action of a suitably 
engineered environment. 
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